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Abstract 

The phenomenon of superscahng for quasielastic lepton induced reactions at energies of a few 

GeV is investigated within the framework of the relativistic impulse approximation. A global 

analysis of quasielastic inclusive electron and charged-current neutrino scattering reactions on 

ffv , nuclei is presented. Scaling and superscaling properties are shown to emerge from both types 

^^ ' of processes. The crucial role played by final state interactions is evaluated by using different 

(N' 

^^ ■ approaches. The asymmetric shape presented by the experimental scaling function, with a long 

^, 

^D . tail in the region of positive values of the scaling variable, is reproduced when the interaction 

\0 ■ 

^^ ' in the final state between the knockout nucleon and the residual nucleus is described within the 
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relativistic mean field approach. The impact of gauge ambiguities and off-shell effects in the scaling 
function is also analyzed. 

PACS numbers: 25.30.Pt; 13.15.+g; 24.10.Jv 

Keywords: Electron scattering reactions, charged-current neutrino reactions, scaling, superscaling, quasielas- 
tic peak, relativistic mean field 



I. INTRODUCTION 

Scaling is a very general phenomenon [l| occurring in various areas of physics that deal 
with probes weakly interacting with many-body systems in which a single constituent in the 
target sj^tem absorbs the eriergy and momentum transfer. The validity of the concepts of 

1 electron scattering at 
. From an exhaustive 



scahng [2] and superscaling [J] applied to inclusive quasielastic fC 

intermediate to high energies has been investigated at depth in J3|, | ^ 

analysis of the (e, e') world data one concludes that the scaling behaviour is highly fulfilled j^. 
One may distinguish between scaling of the first kind, which corresponds to reduced cross 
sections being independent on the momentum transfer g, and scaling of second kind, namely, 
no dependence on the nuclear species. The analysis of data shows that scaling of first kind is 
reasonably well respected at excitation energies below the QE peak, usually called the scaling 
region, whereas scaling of second kind is excellent in the same region. The simultaneous 
occurrence of both kinds of scaling is named superscaling. At energies above the QE peak, 
where nucleon resonances are important, both types of scaling, first and, to a lesser extent, 
second are broken. Scaling violations are shown to reside mostly in the transverse response, 
but not in the longitudinal which appears to superscale J5|- These results are in accordance 
with the important contributions expected in the transverse channel due to effects beyond 
the impulse approximation: inelastic scattering [y, M , corre^ 
currents (MEG) in both the Ip-lh and 2p-2h sectors [d, Ig Ld, lUl IjJ • 

The scaling analysis of QE (e, e') data was extended into the A region |lj|, leading to 
the extraction of two different scaling functions which embody the nuclear dynamics in the 
two regions. The scaling approach has been exploited to predict inclusive charged-current 
(CC) neutrino-nucleus cross sections. This strategy is based on the assumption that a 
universal scaling function exists, which is valid for both electron and neutrino scattering 
reactions, provided that the corresponding kinematics are similar. This hypothesis of a 
universal scaling function, which is true by construction in the relativistic Fermi Gas (RFG) 
model [1^, was further investigated when final state interactions (FSI) were included. The 
analysis performed in |l5| within the context of a non-relativistic mean field calculation, 
but incorporating important aspects of relativity at the level of the current operators, called 
semirelativistic (SR) approach, proved that scaling properties were highly fulfilled for the 
electromagnetic responses at intermediate to high energies. Moreover, the scaling function 



ations and meson exchange 



extracted from these coincides with that from the (z/, /x) reaction. However, it presents 
a symmetric shape which is not supported by the analysis of (e, e') data. This is not 
unexpected because additional dynamical effects, which are beyond the non-relativistic mean 
field picture considered in [13], are needed in order to reproduce the asymmetry extracted 
from the experiment. 

An investigation of the QE scaling properties of CC neutrino-nucleus scattering within 



m- 



the context of the relativistic impulse approximation (RIA) has been presented in 
Although resorting only to one-body excitations, the RIA has been shown to provide the 
required asymmetry of the scaling function when strong relativistic potentials are included 
in the model. This makes an important difference with previous non-relativistic, or SR, 
calculations based on the impulse approximation [l5|, [13, ll8[. The superscaling function 
evaluated from QE {u, fi) calculations was compared with the (e, e') phenomenological one, 
showing the capability of RIA models to yield the required properties of data. 

In this paper we extend the investigation on scaling by performing a global analysis of 



(e, e') and (z/, /i) reactions within the RIA 
scahng and superscaling studies p, llJ, bj 



ramework. We follow the general procedure of 
la|. First, we calculate inclusive cross sections 
within a specific model and then, obtain scaling functions by dividing them by the relevant 
single-nucleon cross sections weighted by the corresponding proton and neutron numbers p, 
|l9|. The scaling function so obtained is plotted against the scaling variable ip{q,uj), and its 
scaling properties analyzed, i.e., we explore its dependence on the transfer momentum (first 
kind scaling) and on the specific target nucleus (second kind scaling). 

A detailed study of the off-shell effects and gauge ambiguities in the scaling function is 
also presented. The analysis performed and its comparison with data give us important 
clues on the validity of the different theoretical descriptions considered. Furthermore, the 
consistency between (e, e') and (z/, /i) calculations, which reflects the universality property of 
the superscaling function, is also clearly illustrated. Finally, we show that the RIA approach, 
in spite of its simplicity, gives rise to the required shape presented by the expermental 
scaling function. Scaling and superscaling ideas have been carried a step further to include 



neutral-current (NC) neutrino-nucleus scattering processes in ^0|. Here, NC differential 
cross sections were obtained by making use of the phenomenological (e, e') scaling function, 
i.e., the universality property was assumed. It will be very interesting to investigate scaling 
for NC reactions within the RIA framework. This will allow us to prove if various RIA-FSI 



models superscale, and moreover, if the universal asymmetric scaling function emerges also 
from the RIA calculations on NC processes. 

The paper is organized as follows. In Section II we present a brief summary of the basic 
formalism involved in describing inclusive QE electron-nucleus and CC neutrino-nucleus 
scattering processes. We restrict ourselves to the plane wave Born approximation (PWBA) 
and assume the RIA. The analysis of scaling and superscaling with the general expressions 
implied is discussed in Section III. Here we also show the experimental scaling function 
together with a phenomenological fit and a comparison with the simple RFG result. In 
Section IV we present our discussion of the results. We start with a global analysis of 
QE (e, e') reactions where off-shell effects in the differential cross sections are investigated 
at depth. We continue with the study of scaling properties showing that our theoretical 
results do scale even when strong relativistic potentials are present. A separate analysis of 
the different channel contributions is also presented. Comparison with data shows that our 
model calculations do agree with experiment for specific descriptions of FSI. To conclude, 
we compare the scaliiig functions evaluated from (e, e') calculations with those obtained 
from (z^, /u) reactions |la|. Results show that they almost coincide, hence the model is 
consistent with the fulfilment of the universality property. Finally, in Section V we present 
our conclusions. 

II. INCLUSIVE QUASIELASTIC LEPTON SCATTERING FORMALISM: THE 
RELATIVISTIC IMPULSE APPROXIMATION 

This work deals with lepton induced reactions at energies of a few GeV and QE kinemat- 
ics. In particular, we focus on inclusive electron scattering and CC neutrino (antineutrino) 
scattering on nuclei and assume the Born approximation (BA). The leptonic variables (in 
the laboratory system) involved in the processes are K^ = {e, k) the 4-momentum of the in- 
cident lepton (e or 1/^) beam, and K''^ = {e', k') the 4-momentum of the scattered lepton (e' 
or yu). The process is mediated by the exchange of a virtual photon (electron scattering) or a 
charged vector boson (CC neutrino scattering) with 4-momentum Q'^ = {K — K')^ = {u, q). 



The general 



works y, S ll 



brmalism for (e, e') and {u, fi) reactions has been presented in previous 



H 



2l|. Here we simply summarize those basic aspects needed for later 



discussion. Assuming PWBA, i.e., one virtual particle exchanged and leptons described as 



free particles, the QE differential cross section can be expressed in terms of separate nuclear 
response functions. In the case of (e, e') reactions we may write 

da 

J (e,e') 

where Q' is the scattered electron solid angle and the term aM represents the Mott cross 
section. Analogously, for CC neutrino scattering reactions the differential cross section can 
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VLR^{q,uj) + VTR^iq,uj] 
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be written in the form 
da 



m 



m 



de'dVt' 
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vccR''^ + 2vclR'^'' + vllR'-'' + vtR' + 2xvt'R 



jT' 



(2) 
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with {e', Q') the muon kinematical variables. The symbol x specifies neutrino-induced re- 
actions (x = +) or antineutrino-induced reactions (x = — ) and the term ctq depends on 
the Fermi constant and the Cabibbo angle (see 13] for its explicit expression). The kine- 
matical factors Vk and Vk come solely from the electromagnetic and weak leptonic tensors, 
respectively, and their explicit expressions can be found in |4Ll8Lll3J. 

The electromagnetic R^ and weak R^ response functions contain the whole dependence 
on the nuclear vertex coupling and are expressed by taking the appropriate components of 
the nuclear tensor j^ |8|, |l3|. This involves the matrix elements of the virtual photon or 
charged boson interaction with the nuclear electromagnetic or weak current. The inclusive 
hadronic electromagnetic tensor reads 



W^'^iq,co) = Y$SiEj -E,- uj){f\Ji:^{Q)WU\J:^{Qm 



(3) 



where \i) describes the initial target state and |/) represents a specific many-body final 
nuclear state. The term JemiQ) refers to the nuclear electromagnetic many-body current 
operator. A similar expression to (jH)) should be written for the weak tensor W^^'^ in terms 
of the nuclear weak many-body current operator J^{Q). The electromagnetic tensor given 
in Q is an exceedingly complicated object which includes all possible final states that can 
be connected with the initial ground state through the action of the many-body current 
operator. 

In this paper we restrict ourselves to the QE kinematic regime and we adopt the rela- 
tivistic impulse approximation. Within the RIA the many-body nuclear current operator is 
simply given as a sum of single-nucleon current operators that only couple the target ground 



state to scattering states lying in the one-body knockout space. The RIA approach 



extensively applied in investigations of exclusive electron scattering reactions 
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las been 



24 



Further details on the model for neutrino-nucleus scattering reactions have been presented 



in mi^ 



271 



28l |. Within the RIA framework the main ingredient needed to evaluate the 



electromagnetic and weak tensor is the single-nucleon current matrix element, 



(>(g))= /rfre^''->^(pp,r)f^^r(r). 



(4) 

Here ?/';g"(r) and ippiPF, r) are the wave functions for the initial (bound) nucleon and for the 
emitted nucleon, respectively, and F^ is the corresponding single-nucleon current operator 
for electron (F^^) or weak CC neutrino (F^^) scattering. 

We describe the bound nucleon states as self-consistent Dirac-Hartree solutions, derived 
within a relativistic mean field (RMF) approach using a Lagrangian containing a, u and p 
mesons J2^ |2l| ■ The outgoing nucleon state is described as a relativistic scattering wave 
function. Different options have been considered. First, the Relativistic Plane- Wave Impulse 
Approximation (RPWIA), namely, the description of the knockout nucleons by means of 
plane-wave spinors. Second, the effects due to FSI between the ejected nucleon and the 
residual nucleus. In our model FSI effects are described by using Dirac equation solutions in 
presence of relativistic potentials. This constitutes the Relativistic Distorted- Wave Impulse 
Approximation (RDWIA) J22|. 

The use of energy-dependent complex relativistic optical potentials fitted to elastic proton 
scattering data has proven to be successful in describing exclusive (e, e'p) scattering reactions 



under QE kinematics 



m, 



H 



E 



25l . l31| . In this case of exclusive reactions, the optical 



potentials are built to reproduce the contribution from the elastic channel. For inclusive 
processes such as (e, e') and (i^, /i), the contribution from the inelastic channels should be 
retained. Ignoring them would lead to an underestimation of the inclusive cross section 21, 



33[. Multiple nucleon knockout effects have been treated in detail within the context of 



the Green function method 



M 



35, 



A simple way of obtaining the inclusive strength 
within the RIA is to use purely real potentials. We consider two choices for the real part. 
The first uses the phenomenological relativistic optical potential from the energy-dependent, 
A-independent parameterizations (EDAIC, EDAIO, EDAICa) derived by Clark et al. J38l |. 
but with their imaginary parts set to zero. The second approach consists of describing the 
outgoing nucleon by means of distorted waves obtained with the same relativistic mean field 



used to describe the initial bound nucleon states. We refer to these two FSI descriptions as 



real Relativistic Optical Potential 



Green function techniques 
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rROP) and RMF, respectively. Dispersion relation and 
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37( 1 lead to results which are close to those obtained 



in the impulse approximation with either the rROP 



or the mean field 
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Concernin g; th e current operator we make use of the relativistic free nucleon expres- 



sions 



Q 
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40|. For electromagnetic (e, e') processes the three usual options, denoted as 



CCl, CC2 and CC3, are considered: 
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(5) 
(6) 
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-ip(n) 



where Ff and F2 are the Pauli and Dirac proton (neutron) form factors, respectively. 



that depend only on Q , and the on-shell 4-momentum P = {E, p) with E = Jp^ + m,% 
and p the bound nucleon momentum has been introduced. Note that the three operators 
are equivalent for free on-shell nucleons (they are connected by the Gordon transformation). 
However, the RIA deals in general with off-shell bound and ejected nucleons. Hence the 
three operators lead to different results. Moreover, the current is not strictly conserved and 
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m 



uncertainties dealing with the election of gauge also occur 

The relativistic charged weak current of the nucleon is given as P^ = Ty — T'^, where the 
vector and axial-vector current operators read 
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Fii^ + ii^^T^'Qu 
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Gai" + ^Q^ 
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(8) 
(9) 



with FY2 the isovector nucleon form factors given in terms of the electromagnetic ones as 



WV _ T?P 

-^1,2 



F12 ~ Pi 2- The axial- vector and pseudoscalar form factors are parameterized as 

9A 



Ga — 

Gp = 

with gA = 1.26 and Ma = 1032 MeV (see 



1 - Q^Ml 



Am% 



mt 



Q- 



:G 



A 



(10) 

(11) 



44 



III. SCALING AND SUPERSCALING AT THE QUASIELASTIC PEAK 

Detailed studies of scaling and superscaling for electron-nucleus cross section have been 
presented in J3|, y, |5[ . The analysis of the (e, e') world data has shown the quality of the 
scaling behavior: scaling of the first kind (no dependence on the momentum transfer) is 
quite good at excitation energies below the QE peak, whereas scaling of second kind (no 
dependence on the nuclear species) works extremely well in the same region. In this paper 
our aim is to investigate the QE scaling properties of electron-nucleus and CC neutrino- 
nucleus scattering within the context of the RIA. Assuming various RIA models we prove 
that they do superscale, and we compare the associated scaling functions with the (e, e') 
phenomenological one. 

In what follows we present the basic expressions needed to get the scaling functions. 
Several choices have been proposed in the literature for the appropriate scaling variable. 
Here, following the analysis of the RFG model, we adopt the dimensionless variable denoted 
as^p'{q,uJ), 

^'^^ , ''''' (12) 

^^ ■ I{1 + X')t' + kJt'(1 + t') 



where X' = {to — Eshift)/'2'm'N, i^ = q/'^ttT'N, t' = k^ — A'^, and ^i? = w 1 + {kp/m^y — 1. 
The term kp is the Fermi momentum and the energy shift Eghift, taken from p|, has been 
introduced to force the maximum of the cross section to occur for ip' = 0. As usual the 
notation ip refers to the scaling variable when E shift = 0. 

For inclusive QE electron scattering processes, the superscaling function is evaluated by 
dividing the differential cross section ((H) by the appropriate single-nucleon eN elastic cross 
section weighted by the corresponding proton and neutron numbers j^, y, ll2| involved in 
the process. We may write 

da 



fii^^q) 



de'dQ' 



ice') 



(13) 



o-M [ViGLiq, uj) + VrGriq, uj)] ' 
The scaling behaviour can be also analyzed by taking into account the separate electromag- 
netic longitudinal (L) and transverse (T) contributions. Thus the following scaling functions 
are introduced: 

8 




-1.5 -1 



-0.5 0.5 1 1.5 2 



FIG. 1: RFG superscaling function compared to data and a parameterization of the results. 



Mi^'.q) = k, 



R^iq.u) 



GTiq,io) 
The single-nucleon functions Gl and Gt are given by 

[nyr) \gI + W2A 



(15) 



Gl = 
Lxj' = 
where the function A reads 

As usual one has 



2tGI, + W2A 
2K[l + ip{l + ,lj^)/2] 



(16) 
(17) 



t(1 + t) 



T 



K 



^ui'r)-j 



K^ 



(18) 



Gp = ZGpp + NGp^ , Gj^ = ZGj^p + NGj^^ 



Wo 



1 + r 



Gjp + tG, 



M 



(19) 



involving the proton and neutron form factors weighted by the proton and neutron numbers 
Z and A^, respectively. 
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At sufficiently high energies the function / depends only on the scaling variable if)' but 
not on the transferred momentum q. Moreover, /(^') becomes also independent on the 
momentum scale in the problem, that is, independent of kp. The scaling behavior has been 
clearly demonstrated from the analysis of the QE (e, e') world data J3l,l4|. The investigation 
of the separate contribution of the longitudinal and transverse response functions has shown 
that scaling violations occur mainly in the T-channel due to significant effects introduced 
by MEC, correlations and inelastic scattering. From the global analysis presented in [3|, y] a 
universal ea;]5enmenta/ superscaling function fexpi'^P') has been defined. In Fig. ^ we present 
fexpi'4'') averaged over the nuclei employed in the analysis, together with the corresponding 
fit. As noted, the experimental scaling function presents an asymmetric shape with a tail 
that extends towards positive values of the scaling variable ip'. This is in contrast with the 
RFG superscaling function given by fRPci^P') = (3/4)(l — ip''^)6{l — ip'^), that is symmetric, 
limited strictly to the region \ip'\ < 1 and with a maximum value of 3/4. 

The general analysis of superscaling for CC neutrino-nucleus scattering has been pre- 
sented in previous works (l3l . ll5L llfil | . Here the superscaling function is obtained by dividing 
the differential cross section evaluated within the RIA by the corresponding weak single- 
nucleon cross section as given explicitly in Eqs. (15,45,52,86-94) of [l3|. Details and specific 



expressions are also given in appendix C of 



15| . This theoretical scaling function calculated 



from CC neutrino-nucleus cross sections can be compared directly with the one correspond- 
ing to (e, e') scattering calculations as well as with fexpi'^P') shown above. This allows a 
check on the universality assumption of f{ip') and on the capabilities of different RIA mod- 
els to yield or not the required properties of the experimental scaling function. By analogy 
to (e,e'), and in addition to the scaling function obtained from the CC neutrino-nucleus 
cross section, one may also construct the separate contributions given by the longitudinal L, 
transverse T and axial- vector transverse T' responses, fL,T,T'{'4'')- Here the L weak response 
includes the contribution from the three terms in (j2)), namely, vccR'"'~^ + '2'<JclR'"^ + 'VllR^^- 
In next section we present a detailed study of scaling and superscaling properties for (e, e') 
and (z/, n) reactions within the scheme of the RIA. 
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IV. RESULTS 



In this work we consider the PWBA, i.e., Coulomb distortion of the leptons is neglected. 



21 



Checks made for light-to-medium nuclei within the effective momentum approxmation 
27l | show that these effects are within a few percent for the high energy lepton kinematics 
considered in this work. In this sense, our general conclusions about scaling are not modified 
by them. Obviously, caution should be drawn on the analysis of heavier nuclear systems 
where a careful description of Coulomb distortion effects in the leptons (electrons and muons) 
involved in both processes (e, e') and {v, fi) is required. 



A. Inclusive QE (e, e') reactions 



1. Differential cross sections 
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FIG. 2: Differential cross section for QE (e,e') on ^■^C. FSI are described within the RMF and 
the rROP models. Results correspond to CCl and CC2 current operators. Also for reference the 
RPWIA result is presented. The incident electron energy is e = 1 GeV and the scattering angle 
^e = 45*^. See labels for reference on the different curves. 

In what follows we present predictions for QE (e, e') reactions on ^^C within the frame- 
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work of the RIA. Results correspond to fixed values of the incident electron beam energy, 
e = 1 GeV and scattering angle, 6^ = 45". The Lorentz gauge has been selected. In next 
section we present results for the scaling function concerning gauge effects. In Fig.|21we show 
the differential cross section evaluated for the two current operators, CCl and CC2, and FSI 
included through the rROP and RMF potentials. As observed, the CCl choice leads to a 
significantly larger cross section, particularly within the RMF approach. Concerning FSI, 
the use of the RMF potential gives rise to a clear asymmetry in the cross section with a 
pronounced tail extending towards higher values of the transfer energy u. For reference we 
also include in Fig. |21the curve corresponding to RPWIA, i.e., no FSI. In such a case, the 
effects introduced by the current operator choice are minor. The difference observed between 
both FSI descriptions at high cj-values is linked to the behaviour of the two relativistic po- 
tentials: whereas the RMF contains strong energy-independent scalar and vector potentials, 
the energy dependence of the rROP makes its scalar and vector terms to be importantly 
reduced for high nucleon kinetic energies (high transfer energy). Hence, rROP results get 
close to the RPWIA ones for large w-values. 

The asymmetry in the differential cross section is proved to be an effect entirely linked 
to the FSI description. Only in presence of relativistic optical potentials with strong scalar 
and vector terms (RMF approach), a significant shift of strength to higher values of u is 
shown to occur. Details on the specific mechanism that produces the asymmetric tail in the 
cross section within the RMF-FSI approach are given in (4(t| . 

The importance of the current operator choice for the two FSI models is further investi- 
gated in Fig. 01 where a separate analysis of the longitudinal and transverse contributions 
to the cross section is presented. As observed, the pure longitudinal response is almost 
identical with both current operators. In the case in which FSI are neglected, this has been 



also found previously |4l| and it was proven as due to the validity of the Gordon transfor- 
mation for the longitudinal contributions. For the case of the transverse channel, the CCl 
contribution is much larger than the CC2 (likewise for CCS) one. Notice that the magnitude 
of this discrepancy depends on the specific FSI description, being larger for RMF, whereas 
in the plane wave limit (RPWIA) both currents lead to very similar results. The ambiguity 
introduced by the current choice for inclusive (e, e') reactions has been already signaled in 
some previous works |47l . |48 1 . 

The large effects introduced by the current operator within the RMF (likewise for rROP) 
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FIG. 3: Longitudinal (L) and transverse (T) contribution to the cross section. Left panel corre- 
sponds to RMF description of FSI and right panel to rROP. In both cases results are presented for 
CCl and CC2 prescriptions. 



approach can be also analyzed by comparing directly theoretical and experimental cross 
sections. This may allow us to determine which particular choice is more appropriate. 
In Fig. m we compare data 



49 



50 



51 



521 I53L I54| with the results corresponding to the 
RIA-RMF approach with CCl and CC2 at very different kinematics. As a general rule 
we conclude that CCl tends to overpredict data whereas the reverse applies to CC2. This 
outcome favours clearly the CC2 option as far as effects beyond the QE peak (A and MEC 
effects) may also play a significant role in the analysis of the data even at the maximum of 
the peak. However, data are not conclusive yet as far as the A and MEC contributions have 
not been evaluated. 



2. Analysis of scaling behavior 

The important effects already shown in the differential cross sections are also visible 
in the scaling function. In the following we rely on the investigation of the superscaling 
properties and present results for the scaling function /(V'') (fT^ . as well as the separate 
longitudinal and transverse contributions, i.e., fii'^'') ^md /t(V'')- ^ comparison between 
the three scaling functions is presented in Fig. El where for simplicity only the RMF-FSI 
model has been considered. Similar conclusions are drawn for the rROP approach. In 
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FIG. 4: Differential cross section for (e,e') reactions on ^^C. Theoretical predictions correspond 
to the RMF-FSI approach with CCl (solid line) and CC2 (dashed line). Different kinematics (see 
labels) have been considered. 

addition to the usual CCl (top panel) and CC2 (middle) prescriptions we also include the 
results for the CCS choice (bottom panel). As observed, the difference between f{ip') and 
the contributions /l(V'') ^ind fri'^P') is very large for CCl, being much smaller for CC2 and 
almost negligible for CCS. This means that zeroth kind scaling is fully broken for CCl and 
RMF, whereas only a mild (negligible) violation is observed for CC2 (CCS). This behaviour 



14 




1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 



FIG. 5: Analysis of zeroth kind scaling. The global scaling function /(V'') (solid line) is com- 
pared with the separate L (dashed) and T (dotted) contributions. All results correspond to RMF 
description of FSI and current operators: CCl (top panel), CC2 (middle) and CCS (bottom). 

of the superscaling function is in accordance with the very diverse contributions given in 
the T channel by the different current operators. In particular, the CCl current and RMF 
description of FSI lead to an important increase of the T channel strength compared with 
the single-nucleon contribution. On the contrary, the longitudinal function fi is found to 
be basically the same for the three choices of the operator (likewise for rROP). 

Scaling of first kind is explored in Fig.lHlwhere we present f{ip') for three different values of 
the incident electron energy, e = 1, 1.5 and 2 GeV. Results are shown for the two different 
descriptions of the FSI: RMF (top panels) and rROP (bottom panels). In each case we 
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FIG. 6: Analysis of first kind scaling. Scaling function for three values of the incident electron 
energy. Kinematics and target as in previous figures. Results correspond to RMF model (top 
panels) and rROP (bottom) and the two choices of current operators: CCl (left) and CC2 (right). 

make predictions for the two usual current prescriptions, CCl (left) and CC2 (right). As 
observed, the scaling function for the rROP model shows a very mild dependence on the 
transfer momentum in both positive and negative ip' regions, i.e., scaling or first kind is well 
satisfied. In the case of the RMF model, a slight shift occurs in the "scaling region" ip' < 0, 
whereas for -i/^'-positive the model breaks scaling at roughly 30%. This violation is not in 
conflict with (e, e') data that indeed leave room for some violation of the first-kind scaling 
in this region, due partly to A production and partly to other contributions, such as MEC 
and correlations. A separate analysis of the L and T channels in the scaling function (see 
Fig. [7j) shows that the breakdown of scaling within the RMF description of FSI is similar 
for both channels independently of which current choice is considered. In the case of rROP 
(and RPWIA) calculations scaling of first kind is excellent for all current operators and in 
both channels. 
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FIG. 7: As in previous figure but for the separate L (top panel) and T (bottom) contributions. Only 
the RMF approach is considered. In the case of the T response, the larger functions correspond 
to the CCl operator and the smaller ones to CC2. In the L channel there is almost no distinction 
between both operators. 



Scaling of second kind, i.e., independence on the specific nuclear system, is analyzed in 
Fig.lSl Here we present tlie results corresponding to CCl (left panels) and CC2 (right panels) 
operators and the two descriptions of FSI: RMF (top panels) and rROP (bottom panels). 
In each case we compare the superscaling function evaluated for three different nuclei: ^^C 
(solid line), ^^O (dashed) and "^"Ca (short-dashed). The values of the Fermi momentum 
considered [sj correspond to kp = 216 MeV/c (^^O), kp = 228 MeV/c {^^C) and kp = 241 
MeV/c (^°Ca). As shown, the effects introduced by changing nucleus are very small for 
the CC2 current choice and the two FSI descriptions. This is in complete accordance with 
data which show that second-kind scaling is excellent. On the contrary, the CCl operator 
leads to a significant breakdown of scaling behaviour which affects both FSI descriptions. 
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FIG. 8: Analysis of second kind scaling. Scaling function for three nuclei: ^^C, ^^O and ^°Ca. 
Kinematics as in previous figures. Results for the RMF (top panels) and rROP (bottom) models 
are presented for CCl (left) and CC2 (right) prescriptions. 

This violation of the second kind scaling conies totally from the transverse response. This is 
clearly illustrated in Fig. 121 where we present results for the separate functions /l(V'') (top 
panel) and fri'^'') (bottom) corresponding to the RMF model. Similar results are obtained 
for the rROP approach. From these theoretical results and the exhaustive analysis of the 
(e, e') world data, which proves the excellent quality of second-kind scaling behaviour, one 
may question the validity of the CCl operator used in describing QE (e, e') processes. 

3. Comparison with experiment 

A comparison between the theoretical superscaling functions and the averaged QE phe- 
nomenological function obtained from the analysis of (e, e') data is presented in Fig.^J We 
have selected the CC2 operator and show results for the three different descriptions of the 
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FIG. 9: Separate L (top) and T (bottom) contributions to second-kind scaling analysis. Results 
correspond only to the RMF model. 



continuum final state, namely, the plane wave limit (dashed line), rROP (dot-dashed) and 
RMF (solid). Lorentz gauge has been assumed. The symmetric character of the RPWIA 
and rROP curves differs clearly from the experimental analysis. On the contrary, the RMF 
approach displays an asymmetric shape with a long tail extended to positive values of the 
scaling variable ip' which follows closely the behaviour of the phenomenological function. As 
already mentioned in previous work [l^ , this asymmetry of the RMF description constitutes 
a basic difference from other models presented in the literature 55|, |56|, |5|^ where the long 
tail in f{'4)') is largely absent. 

The asymmetry in data has usually been ascribed to ingredients beyond the mean field, 
such as short-range correlations, induced non-localities and two-body currents. Within a 
non.eMiv.«c ap.oach .c. .n,.dient= a. needed . o.de. .o .e. ..e ..„.e.y UU 
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FIG. 10: Scaling function for ^^C{e,e') evaluated with the RPWIA (dashed line), rROP (dot- 
dashed) and RMF (solid) approaches compared to the experimental function together with a phe- 
nomenological parameterization (dotted). 



58l |. However, here we show that a large amount of the asymmetry is indeed obtained within 
the framework of the RIA and a RMF description of the final continuum nucleon states. 
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62| where a comparison between 



This is in accordance with some previous works |59l 
Dirac-Brueckner-Hartree-Fock (DBHF) calculations and Dirac-Hartree ones indicates that 
effects from correlations and Fock terms in the DBHF calculation can be accounted for 
by the simple Dirac-Hartree approach fitted to saturation properties of nuclear matter. 
This is at variance with the non-relativistic mean field case. In this respect, note that 
the Dirac equation in presence of scalar and vector local potentials can be reduced to a 
non-relativistic Schrodinger-like equation with energy-dependent and non-local terms 23|. 
Results in Fig. ITUl show that the asymmetry in the scaling function can be produced via local, 
energy-independent relativistic potentials within the impulse approximation, and moreover, 
such asymmetry is very close to the experiment. This outcome does not contradict the 
additional role that may play correlations and exchange currents not accounted for within 
the relativistic mean field calculation. However, the small magnitude of the local central 
and spin-orbit potentials involved in the non-relativistic approach cannot yield a significant 
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FIG. 11: Scaling functions evaluated for different currents and gauges compared to experiment. 
Left panels show /(V'') as well as the contributions /l(V'') a-iid /rlV'O for the three choices of the 
current: CCl (top), CC2 (middle) and CCS (bottom) and the gauge fixed to Lorentz. Right panels 
show the results for the longitudinal scaling function /l('0') for the three currents and the three 
gauges: Lorentz (solid line), Coulomb (short-dashed) and Weyl (dashed). 



asymmetry. Only a strong non-locality of the potentials (effective values of the mass and 
energy) may give rise to an asymmetric differential cross section J46| . 

To complete the analysis, in Fig. ^2 (left panels) we select the RMF description of FSI 
and compare data and the fit curve with the theoretical results for the three choices of 
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the current operator: CCl (top panel), CC2 (middle) and CCS (bottom). In each case 
we also present a separate analysis of the L and T channels involved in the process. As 
already shown in previous results, the longitudinal contribution /l('?A') does not depend 
on the current operator choice and agrees nicely with experiment. On the contrary, the 
strength of the transverse response with the CCl current leads to a function fri'ip') that is 
roughly twice the data. Discrepancies between /l(V'') and fri'ip') are mild (negligible) for 
CC2 (CCS) leading to a global scaling function /{ip') which is in accord in both cases with 
the experiment. All the curves presented for the longitudinal contribution /l(V'') satisfy the 
Coulomb sum rule, i.e., they integrate to unity. 

Up to present we have only discussed calculations corresponding to the Lorentz gauge. 
Results for the Coulomb gauge are similar for all FSI descriptions whereas the Weyl gauge 
leads in most of the cases to very important discrepancies. This is clearly illustrated in the 
right panels of Fig. ITT] where we show the longitudinal scaling function /l(V'') evaluated 
for the three gauges: Lorentz (solid line). Coulomb (short-dashed) and Weyl (dashed), and 
the three current operator choices: CCl (top panel), CC2 (middle) and CCS (bottom). 
The RMF approach for the final state has been assumed. The discussion of results follows 
in general similar trends for the rROP model. Note that the gauge does only affect the 
longitudinal response. As observed in Fig. [TTl Lorentz and Coulomb gauges lead to very 
close results for all currents, particularly for CC2, where there is no distinction between 
the two curves. In the case of the Weyl gauge, an important difference emerges between 
the CC2 operator and the other two, CCl and CCS. In the former, the scaling function 
/l(V'') is almost the same for the three gauges. This is in accordance with the fact that the 
continuity equation is highly fulfilled within the RMF model and the CC2 current operator. 
By contrast, the CCl and CCS operators within the Weyl gauge give rise to longitudinal 
scaling functions which depart significantly from data being much larger (smaller) for CCl 
(CCS). In both cases the Coulomb sum rule is clearly violated. These results reinforce our 
confidence in the adequacy of descriptions of inclusive (e, e') reactions when based on the 
RMF-FSI approach and the CC2 current operator. Only in this case the results do not 
depend on the specific gauge selected^ and moreover, they are also shown not to be modified 



^ Although not shown, the Weyl gauge leads to different results when the rROP model is assumed indepen- 
dently of the current operator selected. 
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B. Inclusive QE charged-current (i/, /x) reactions 
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FIG. 12: Consistency of scaling functions evaluated for (e, e') and (i^, /x) reactions. The incident 
electron (neutrino) energy is fixed to 1 GeV and the scattering angle to 45'^. Left panel refers to 
results obtained with the CC2 current (electron scattering) and a separate L and T contribution 
is presented. Right panel shows results for CCS choice. 



The analysis of scaling and superscaling for CC neutrino-nucleus reactions has been pre- 
sented in 13|, |l6[ . It is important to point out that any reliable calculation on neutrino- 
nucleus cross sections requires first to be tested against electron scattering data. Hence two 
different approaches can be pursued. First, using the scahng behaviour of (e, e') cross sec- 
tions and the universality property of the superscaling function, we can make predictions for 
inclusive (z/, yu) reactions by taking the empirical electron scattering scaling function fexpii^')- 
This strategy, applied not only to the QE regime but also to A kinematical region, was an- 

I i n 

alyzed at depth in jl3]. The second approach, considered in jl6|, consists of evaluating 
explicitly f{ip') for (z/, /x) reactions within a specific model, namely, RDWIA. The scaling 
function obtained in this way can be compared directly with the model predictions given 
for (e, e') processes. This allows us to check the scaling behaviour of the calculations and 
moreover, the consistency of the universality assumption of f{ip') and also the capability of 
the model to reproduce the experimental data. 
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In jl6| we presented a detailed investigation of CC neutrino-nucleus scattering reactions 
within the RIA framework. We proved that superscaling is verified to high accuracy by 
the model calculations even in presence of strong relativistic potentials. Importantly, the 
results obtained when FSI were described by means of the RMF potential presented the 
right asymmetry compared with data. This is fully consistent with the discussion outlined 
in previous section concerning the study of inclusive (e, e') reactions. This consistency is 
clearly illustrated in Fig. ^J where we compare the scaling functions evaluated from (e, e') 
and {u, /i) reactions. We only consider the RMF-FSI case as this is the only model which is 
in accordance with data. However, the consistency between electron and neutrino scattering 
calculations applies also to rROP and RPWIA approaches. In the left panel of Fig. ^J the 
CC2 prescription has been assumed for (e, e') and the separate contribution of both channels, 
L and T, is also shown. For completeness, the case of the CCS operator is considered in the 
right panel where the fit curve to experiment is also drawn for reference. From these results, 
it is clear that the universality assumption of the scaling function is highly fulfilled within 
the present RIA model. This supports the use of the experimental (e, e') scaling function in 
order to predict reliable neutrino- nucleus scattering cross sections 13l |. 

Although not presented for simplicity, a separate analysis of the various channels, L, T 
and T', contributing to (z/, /i) reactions shows that frijp') = fr'ii'') = fi'ip'), i-e., scaling of 
zeroth kind is verified. The longitudinal contribution leads to a scaling function /l(V'') which 
departs significantly from f{ip'). However, one should be cautious because the L contribution 
to inclusive {u, fj,) cross sections is almost neghgible compared with the transverse, T, T', 
ones. 

To conclude with the analysis of results, we present again in Fig. ^J the scaling function 
predictions for both (e, e') and (z/, /i) reactions with the RMF description in the final state, 
and compare with data and the fit curve. Here we use a logarithmic scale in f{ip') in order to 
enlarge the discrepancies between theory and experiment in the scaling region, i.e., negative 
'^/''-values. In this region, our model predictions for electron and neutrino processes are in 
fully agreement and tend to underpredict the data. This is not unexpected because the model 
is entirely based on one-body phase space. Ingredients beyond the impulse approximation, 
i.e., multinucleon knockout, either induced by exchange currents, correlations or rescattering 
effects, are surely needed to get more strength in the scaling function which will be then 
closer to the experiment. This is in fact the case of the Coherent Density Fluctuation Model 
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FIG. 13: As in previous figure but in logarithmic scale. We also include data and compare with 
the predictions given by the CDFM model (see text for details). 



in y, 



(CDFM) for correlations presented in J55|, |56|, |57|. This model is an extension of the RFG 
applied to finite nuclei and its prediction is also shown for comparison in Fig. ^1 The 
CDFM result, compared with the RMF calculations, presents more strength in the negative 
^/''-region being closer to data. However, CDFM is manifestly symmetrical around the QE 
peak.^ As discussed in previous sections, the asymmetry of the scaling function comes 
mainly from the inclusion of FSI in the reaction mechanism. This ingredient, when based 
on the RMF description, gives rise to the right asymmetry shown by the experiment. 

A detailed analysis of scaling and cross sections within a non-relativistic mean field ap- 
proach has been presented in |l^. In Fig. El we compare CC neutrino-nucleus cross sections 
obtained with the RIA model and the semirelativistic (SR) one developed in [15[. The SR 
approach includes important relativistic ingredients in the current operator. This has been 



^ The author is aware of a new development of the CDFM model in which asymmetry is incorporated in 
an effective way \6^ . 
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FIG. 14: Predictions for ^'^C{v,n) differential cross sections. Tlie RIA model is compared with the 
SR approach. Two different descriptions of FSl have been considered. First the plane wave limit 
which leads to very similar results. Second FSI described by means of the RMF model (dashed) 
and making use of a Woods Saxon potential in the SR framework (dotted). 



already illustrated in previous works J9|, |l5|, |6J, |65[ and it is also clearly shown up by results 



in Fig. El corresponding to the plane wave limit; the fully relativistic and the SR calcula- 
tions lead to very similar curves in spite of using different descriptions of the initial bound 
nucleon wave functions: solutions of the Dirac-Hartree model and the Schrodinger equation 
with a Woods Saxon potential, respectively. On the contrary, relativistic and SR results are 
very different when FSI are included by using potentials (RMF and Woods Saxon) in the 
final state. Whereas the RMF shows a significant asymmetry, the SR-FSI approach shifts 
the SR-PWIA result to higher muon energies, but maintains the global symmetry of the 
cross section p, ll^ • Asymmetry within the non-relativistic framework requires ingredients 
in the final state beyond the impulse approximation and mean field picture: correlations, 
non- localities, exchange terms, multinucleon emission. By contrast, most of the asymmetry 
involved by the experiment is already obtained within a simple, fully relativistic mean field 
model based on the impulse approximation and describing FSI by means of strong local 
scalar and vector potentials. This makes a crucial difference between both, relativistic and 
non-relativistic, descriptions. 
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V. CONCLUSIONS 

In this paper we have presented a global study of scaling and superscaling properties 
concerning inclusive QE electron and CC neutrino-nucleus scattering reactions. The general 
framework in which the calculations have been performed is the RIA. This is a simplified 
description of the reaction mechanisms because it presupposes that the processes to be 
calculated depend entirely on one-body phase space. However, in spite of its simplicity, the 
RIA has been used to describe successfully QE (e, e'N) reactions. For the inclusive (e, e') and 
{v, fi) processes we are interested in, a crucial ingredient is the description of FSI. We have 
considered different options consistent with retaining in the calculations the contribution 
from the inelastic channels. This is provided in our model by using purely real potentials. 
Three different cases have been explored: i) the plane wave limit, ii) using the real part of 
the energy- dependent relativistic optical potentials parameterized by Clark et al. [38|, and 
iii) employing the same relativistic mean field potential considered in the description of the 
initial bound nucleon states. 

The differential (e, e') cross sections are shown to be very sensitive to some ingredients 
of the model, particularly, the specific current operator selected (CCl vs CC2) and the 
relativistic FSI description. The former, off-shell effects, enter essentially in the transverse 
channel, whereas the longitudinal one shows a very mild dependence with the current op- 
erator choice. Concerning FSI, a basic difference emerges when comparing the rROP and 
RMF approaches to the cross section. The latter presents a significant asymmetry extended 
to higher values of the transfer energy. The orig in of this effect is directly linked to the use 
of very strong relativistic optical potentials 46| . 

The main goal of this paper is centered in the analysis of the scaling behaviour of the cross 
sections and separate responses. Our results show that scaling of zeroth order is verified for 
the CC2 and CCS current operators; however, the CCl leads to an important breakdown 
of scaling because of the extremely large contribution given by the transverse response. 
This outcome is strictly true when FSI effects are taken into account. In the plane wave 
limit (RPWIA), results present a mild dependence with the current operator. Concerning 
first kind scaling, i.e., independence on the transfer momentum, all our results satisfy this 
property which is in accord with the general behaviour of data. Being more precise, only 
in the case of FSI described with the RMF potential, the superscaling function shows some 
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dependence with the transfer niomentum. However, this scahng breakdown is compatible 
with (e, e') data. A separate analysis of the two channels involved in the calculations leads 
also to similar conclusions. 

A special comment should be made concerning second kind scaling, i.e., independence of 
the nuclear species. Recent analyses of (e, e') data conclude that this property is excellent, 
showing only a small scaling violation in the region above the QE peak (positive ip'-vaAnes). 
This means that any theoretical model which does not satisfy this behaviour should be 
dismissed, or at least, clearly disfavoured. In our model, the analysis performed shows that 
the use of the CCl operator leads to a visible breakdown of second kind scaling. This applies 
to both FSI descriptions, rROP and RMF, and it comes totally from the transverse response. 
These results strongly favour the use of the CC2 (likewise CCS) operator, which agrees with 
some general comments made in the past_ concerning the analysis of exclusive QE (e, e'p) 



24| . Note that second kind scaling is highly fulfilled 



data at high missing momentum values 
for the longitudinal contribution fiijp')- 

A comparison with the phenomenological scaling function extracted from (e, e') data 
shows that only the FSI description based on the RMF approach leads to the right amount 
of asymmetry. Within the one-body context, this significant shift of strength to positive 
'^/''-values only occurs in presence of strong scalar and vector relativistic optical potentials. 
This is the case of the energy-independent RMF approach. On the contrary, the energy- 
dependent scalar and vector terms in the rROP are importantly reduced for increasing 
outgoing nucleon energies (positive -i/^'-region) . Comparison with data indicates that the 
RMF-FSI model incorporates important dynamical effects missed by other models. In a 
non-relativistic approach, the asymmetry can be only generated by including ingredients 
beyond the impulse approximation. Here we show that a simple local, energy-independent 
relativistic mean field potential takes care of the basic behaviour presented by (e, e') data. 

The gauge analysis performed in this work leads to an important conclusion concerning 
the RMF-FSI approach. Notice that gauge ambiguity only affects the longitudinal contri- 
bution fii'ip')- Our study shows that Coulomb and Lorentz gauges give rise to very similar 
results, being in general significantly different from the Weyl ones. However, in the particu- 
lar case of FSI described with the RMF potential and the CC2 current operator selected, the 
longitudinal scaling function does not show sensitivity to the gauge election. Only in this 
situation the results of J'l are very similar for the three gauges and, more importantly, all of 
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them essentially satisfy the Coulomb sum rule, i.e., they integrate to unity. By contrast, for 
CCl and CCS prescriptions, the Weyl fi results deviate significantly from the other gauges 
breaking also the Coulomb sum rule. This general result gives us confidence in the adequacy 
of descriptions of QE (e, e') reactions when based on the RMF-FSl approach and the CC2 
current operator. This choice, in addition to be mostly insensitive to gauge ambiguities and 
dynamical relativistic effects, leads to the correct asymmetry shown by data. 

To conclude, we have compared the scaling function evaluated from (e, e') calculations 
with that obtained from (z/, /x) reactions. The results for both processes are very close, and 
this is fully consistent with the general statement on the universality property of the scaling 
function. Within the RIA framework, the superscaling function is basically the same for the 
two t-channel processes considered. This supports the general approach considered in (l3|, 
i.e., the use of the scaling function extracted from (e, e') data to predict neutrino-nucleus 
cross sections. It will be very interesting in the future to analyze if the universality prop erty 
emerges also from RIA calculations applied to NC neutrino-nucleus scattering reactions J2i| . 
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